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1 Overview 

Getting rehable estimates of coronal metalhcity (Z) from X-ray spectra 
obtained with instruments such as ASC A/SIS and Chandra/ ACIS is very 
difficult, because the sole determinant of Z is the ratio of line to continuum 
fluxes, which is not well-determined for low-resolution spectra. Here we 
propose new Bayesian methods which directly model the Poisson nature of 
the data. Our model also accounts for the Poisson nature of background 
contamination, blurring due to instrument response, and the absorption 
of photons in space. The resulting highly structured hierarchical model is 
fit using the Gibbs sampler, data augmentation, and Metropolis-Hasting. 
We demonstrate our methods with the X-ray spectral analysis of several 
apparently coronal metal abundance deficient (" MAD" ) stars. 



2 A Poisson Spectral Model 

The model is designed to summarize the relative frequency of the energy 
of photons arriving at a detector. We model the photon counts in each 
bin as independent Poisson random variables. Specifically, we model the 
high energy tail of the ASCA spectrum (2.5-7.5 keV) as a combination 
of a Bremsstrahlung continuum and ten narrow emission lines, included 
at positions of known strong lines. This source model is combined with 
instrument response, the effective area of the instrument, and background 



contamination to model actual observed counts (van Dyk et ai, 2001). 

Statistical analysis is based on two observations of the same source and 
two of the same background. We use sequential Bayesian analysis for the 
two source observations; the posterior distribution from the first analysis 
is used to construct an informative prior for the second. Non-informative 
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priors were used in the first stage. Model fitting proceeds by using the 
EM algorithm to check for multimodality in the posterior. A Markov chain 
Monte Carlo algorithm is constructed to sample from that posterior. Three 
chains are used at each step of the analysis in order to assess convergence. 
Posterior inference is based on the second half of the draws of all the chains. 
The sensitivity of our results to the choice of prior was investigated by 
altering the prior. We used residual plots for both source and background 
to diagnose the fit of our models. 



3 Results 

We have measured the coronal metallicity for 4 stars, a Aur (Capella), a 
Gem, YY Gem, and II Peg. Our result for Capella [Z = 0.73^q'24) are 
consistent with the independently determined value oi Z = 0.57 — 0.78 



(Brickhouse et al, 200C). Based on our analysis of a Gem {Z — O.SI^q gg 



we find that contrary to classical analyses (Z — 0.25) it is not a MAD star 



(photospheric Z = 0.6, Randich et al. (1994)). YY Gem is found to have 
sub-Solar metall icity (Z = 0A6<^^ ll) b ut higher than EUVE/SW measure- 



ments {Z - 0.1; [Kashyap et al. (1998| )). For II Peg, we find Z = 1.1 



<1.3 
>0.88' 

higher by a factor ~ 10 than preliminary results from Chandra/HETG; this 
discrepancy may be due to anomalies in the abundances of other elements, 
or to a high-temperature component to the emission measure. 
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Abstract 

Metallicity can be measured by analyzing the spectra in the X-ray region and 
comparing the flux in spectral lines to the flux in the underlying Bremsstrahlung 
continuum. In this paper we both explore the sampling properties of the standard 
astrophysical estimates in the small-sample regime appropriate for the data sets of 
interest and propose new Bayesian methods which directly model the Poisson nature 
of the data and thus are expected to exhibit improved sampling properties. Our 
model also accounts for the Poisson nature of background contamination of the 
observations, image blurring due to instrument response, and the absorption of 
photons in space. The resulting highly structured hierarchical model is flt using the 
Gibbs sampler, data augmentation and Metropolis-Hasting. We demonstrate our 
methods with the X-ray spectral analysis of several "Metal Abundance Deflcient" 
stars. The model is designed to summarize the relative frequency of the energy of 
photons (X-ray or 7-ray) arriving at a detector. Independent Poisson distributions 
are more appropriate to model the counts than the commonly used normal 
approximation. We model the high energy tail of the ASCA spectrum of each of the 
stars as a combination of a Bremsstrahlung continuum and ten narrow emission 
lines. The ten emission lines are included at positions of known strong lines. 
Statistical analysis is based on two source observations and one background 
observation. We use sequential Bayesian analysis for the two source observations; the 
posterior distribution from the flrst analysis is used to construct a prior for the 
second. Sensitivity of the flnal results to the choise of prior is investigated by 
altering the prior. 



A Bayesian Poisson Spectral Model 



We use van Dyk et al.'s (ApJ, 2001) hierarchical model for spectral 
analysis of high resolution X-ray and 7-ray spectra. 

Bayesian methodology enables direct modelling of the Poisson nature of 
the source counts and background contamination, instrument response, 
absorbtion and varying effective area. 

We use a generalized linear model for the source continuum (e.g., 
Bremsstrahlung continuum) and allow for several emission lines. 
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Model Specification 



1. Independent Poisson distributions are used to model the 
unobserved uncontaminaited counts. 

Yj ~ Poisson [Xjie)] 

We parameterize the intensity at energy Ej, as a mixture of a 
continuum term and K Gaussian distributions, 

K 
k=l 

f{9^, Ej) represents the continuum and is typically a log linear 
model. 

2. Because of degradation due to the background, the detector, and 
absorption, we model the observed counts as independent Poisson 
variables 

y/''^ - Poisson [^i{e)] 

with parameters 

J 

m = Y.Mij^ji^)djgid^,Ej)^Xfie% 



source 
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where 

• The first term represents the source and the second term 
represents the background. 

• Mij is instrument response, the probability that a photon with 
energy Ej is being detected in observed bin I. 

• dj is the effective area of bin j, normalized to sum to 1. 

• g{9^,Ej) is the probability a photon is not absorbed. 

• \f{9^) is the expected counts due to background. 



4 



Spectra^^mlysi^ofJ 



Metal Abundance Deficient Stars 



We model the high energy tail of the ASCA spectrum (2.5 - 7.5 kev) of 
each of several stars as a combination of a Bremsstrahlung continuum 
and ten narrow emission lines. 

1. The continuum is parameterized as 

<./^7 ^ N Norm TP ih ^ 
f{e\Ej) = — e-^^'^^^ 



2. The ten emission lines are included at positions of known strong 
lines [fi^ are fixed, and their widths are fixed at zero). 

3. The background counts are modeled as: 

- Poisson(6lf ) 

where the background intensities Of ^ are modeled using a second 
log linear model, in particular a Power law, 

\ogef = a-(3\ogEi. 
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The Stars \ 


We analyze the ASCA data of four stars: Capella, a Gem, YY Gem, 


and II Peg. 




The motivation for choosing these specific stars is: 


• To validate our method 


comparing our results with known facts 


about certain stars (Capella, Il-Peg). 


• To compare our method to previous ones and point out the 


modeling advantages {a- 


Gem). We are also planning to analyze a 


few other stars such as ^ Uma and ^ Boo. 


• To get a better understanding of uncommon findings about certain 


stars (YY-Gem). 
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Prior Distributions 



1. We used highly non-informative prior distributions in the first stage 
of the analysis: 



p(log(^)) Pi^ piX^ 

Capella flat flat 7(/x = 1500, a = 1875e3) 

a Gem flat N{n = 5.1, a = 2.5) j{fi = 1500, a = 1875e3) 

YY Gem flat N{fi = 5.1, a = 2.5) 7(/i = 1500, a = 1875e3) 

II Peg flat flat ^{fj, = 150, a = 18750) 



2. The prior for the total count from all the lines, = X^^^i was 
truncated to be less than 70% of the total counts in the selected 
energy range of 2.5 - 7.5 kev; which translates into the following 
truncation points: 



data set 


Capella 


a Gem 


YY Gem 


II Peg 


gis2 


N/A 


2071 


516 


N/A 


sisO 


349 


2470 


627 


2331 


sisl 


284 


N/A 


N/A 


1826 



3. Flat priors were assumed for background parameters. 
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Model Fitting 



Statistical analysis is based on two observations of the same source and 
two of the same background. The background observations are strongly 
informative for 9^ . 

We use sequential Bayesian analysis for the two source observations; the 
posterior distribution from the first analysis is used to construct an 
informative prior for the second. 
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A Sample Statistical Analysis 



The analysis was done by: 

1. Fitting the model using the Markov Chain Monte Carlo algorithm 
(Gibbs Sampler) to sample from the posterior distributions. 

2. Posterior inference is based on the second half of 4000 iterations of 
three Markov chains for a total of 6000 draws for each source 
observation. Each chain uses different, over-dispersed, starting 
values. 

3. Convergence of the MCMC algorithm is verified using the method 
of Gelman and Rubin (1992). 

4. Sensitivity of the results to the choice of prior distribution is 
examined by altering the prior. 

5. The standard deviation of the posterior distribution resulting from 
the first source observation is infiated by a factor of 2 to form the 
prior distribution for analysis in the second source observation. 



9 




a Gem 



Analysis I: 
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-10.0 
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0.5 



0.7 



22 26 30 34 



0.40 



0.55 



95% CI: { 5.43e-005 , 0.000127] 
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Analysis II: 



ln(norm/sqrt(T)) 




o 

d I — . — . — ' — ' — . 

-10.0 -9.4 

95% CI: ( 5.6e-005 , 9.71e-005) 
mean on original scale: 7.39e-00; 



1/(k*Temp) 




0.65 0.80 

95% CI: (1.18, 1.39) 
mean on original scale: 1 .28 



sqrt(Lambda) 




22 26 30 34 

95% CI: (579, 1010) 
mean on onginal scale: 785 



sqrt(omega) 




0.30 0.40 



95% CI: (0.1 , 0.187) 
mean on original scale: 0.142 



Figure 1: Posterior distributions of some of the parameters from the two instruments. 
Lambda is the expected model counts from the hues, and omega is the ratio of the 
total counts in the lines to the total counts in the spectrum. Transformations of 
the parameters are used to produce the distribution nearest to the Gaussian. The 
dashed lines in Analysis I represent the non-informative priors, and the solid lines the 
posterior distribution after having inflated its standard deviation by a factor of two 
to construct the priors for Analysis II (solid lines). 
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Source Model 
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Figure 2: Overview of the fitted model. The first plot compares the fitted model with 
the observed PHA counts. The second one gives the residual for each PHA channel 
after subtracting the background and standardizing by the model standard deviation. 
The third illustrates the stability of the continuum. 
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Figure 3: Those are diagnostics plots for our fitted model. The upper plots are the 
fit of the source model and the corresponding residual plot, and the lower two plots are 
the fit of the background model and its residual plot. The solid lines in the residual 
plots represent plus and minus two standard deviations. 
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Results 



Star 


Data 


Parameters 


Fe/H 






Norm 


kT 


Exp Line Cnt 


%-From Lines 




Capella 


sisl 


0.00442 


0.78 


212 

( 1 97 9Rn^ 


0.26 

(c\ ^ 9/1 n A'\i\ 


-0.137 

(-0.6135, 0.375) 


sisO 


0.00798 

fO Icx ^ 9 9o 9^ 


0.726 

(^u.D, u.yo ) 


177 


0.147 

fn r\fi n 97^ 


a Gem 


gis2 


8.5e-5 


1.52 

/I OCT 1 '70\ 


929 

t^oyu, 1 lou j 


0.26 


-0.08956 

(-0.227, 0.0394) 


sisO 


7.4e-5 

(5.63-5, 9.7e-5) 


1.28 

(1.18, 1.39) 


785 

(579, 1010) 


0.14 

(0.1, 0.19) 


YY Gem 
(3-7.5keV) 


gis2 


0.000613 

(2.4e-4, 1.57e-3) 


1.67 

(1.24, 2.5) 


144 

(91, 204) 


0.178 

(0.098, 0.273) 


-0.3363 

(-0.634, -0.151) 


sisO 


0.00102 

(5.3e-4, 1.9e-3) 


1.44 

(1.18, 1.84) 


101 

(49, 186) 


0.099 

(0.048, 0.174) 


II Peg 


sisl 


0.0037 

(0.00255, 0.0052) 


1.32 

(1.19, 1.48) 


540 

(355, 752) 


0.112 

(0.071, 0.163) 


0.0299 

(-0.056, 0.1135) 


sisO 


0.0028 

(0.0022, 0.0036) 


1.42 

(1.31, 1.54) 


1040 

(872, 1230) 


0.167 

(0.135, 0.199) 



Conclusion 



• Capella coronal Z(Fe) is 0.57-0.78 (Brickhouse et al., 2000, ApJ 539, 
387), which is consistent with our measurements. Our analysis is 
also consistent with Chandra grating data of Capella. 

• (7 Gem XSPEC measurement is [Fe/H] = -0.6, compared to the 
photospheric abundance of -0.3 (Randich, Giampapa, & Pallavicini, 
1994, A&A 283, 893), but we measure the abundance to be 
approximately Solar, and hence conclude that this should not be 
considered a MAD (Metal Abundance Deficient) star." 

• Based on line-to-continuum measurements using EUVE/SW data, 
Kashyap et al. 1998 (AAS 192, 82.01) determined the metallicity of 
YY Gem to be very low, [Fe/H] ~ — 1. This low value may have 
been due to an underestimation of line fiux contribution and an 
overestimation of the continuum level. Our analysis instead shows 
that [Fe/H] is larger, though still sub-Solar. 

• Metallicity derived for H Peg is approximately Solar, which is 
inconsistent with preliminary results from Chandra grating data 
(Drake, J.J., 2001, private communication), which suggest 
metallicities that are almost a factor of 10 lower. The discrepancy 
may be due to anomalies in the abundances of other elements such 
as Si and Mg, or due to a high-temperature component to the 
emission measure. 

• Our findings are inconsistent with 2 or 3 temperature, global fits to 
the spectra (usually carried out over a wider energy range of 0.5-7 
keV). This inconsistency may be attributed to atomic line 
incompleteness at low energies, uncertainty in element abundances. 
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and lack of weight given to the high-energy tail in standard analyses 
becausse of the low counts at these energies. 
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